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Abstract 

The potential between infinitely heavy quarks in a color singlet state is of 
fundamental importance in QCD. While the confining long distance part is in- 
herently non-perturbative, the short-distance (Coulomb-like) regime is accessible 
through perturbative means. In this paper we present new results of the short 
distance potential in coordinate space with quark masses through two loops. 
The results are given in explicit form based on reconstructed solutions in mo- 
mentum space in the Euclidean regime. Thus, a comparison with lattice results 
in the overlap region between the perturbative and non-perturbative regime is 
now possible with massive quarks. We also discuss the definition of the strong 
coupling based on the force between the static sources. 
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1 Introduction 



The potential between two (infinitely) heavy color charges in a singlet state has been 
subject to theoretical investigations for more than twenty years [I], [| ||, In 
the non-perturbative regime it is expected to play a key role in the understanding of 
quark confinement and it is a major ingredient in the description of non-relativistically 
bound systems like quarkonia. In addition it is the basis for the definition of the lattice 
coupling as the potential is given by the vacuum expectation value of the Wilson loop. 

In the perturbative regime it can be utilized to define a physically motivated strong 
coupling which automatically possesses welcome properties such as gauge invariance 
and decoupling of heavy flavors ||. The potential can also be employed for a definition 
of the coupling from the force between the static sources |7j. The latter definition has 
the advantage that the /^-function is unique in that no sign-change of the coupling 
occurs when entering the confinement regime 0. Moreover, the heavy quark system 
is ideally suited to study the overlap region between the non-perturbative and pertur- 
bative treatments of QCD. This latter point is difficult to implement with massless 
dynamical fermions and can only be achieved in the coordinate representation. 

A possibly very interesting application of the two loop mass corrections to the 
heavy quark potential is the effect of a massive charm loop in the MS-bottom mass 
determination ||. Using the physical T-meson for this purpose, the effect of the mass 
shift Srrib depends on (<f) ls \Vp(r,m)\<pi s ), where 4>u denotes the Is ground state wave 
function of the T-meson and Vp the massive fermionic corrections to the potential. 
The effect could be significant if the renormalization scale depends parametrically on 
m c , as indicated by a BLM-analysis ||. For a practical evaluation of these mass shifts 
one needs to have manageable expressions up to the required order in perturbation 
theory. The exact results of Ref. || are not suitable due to the remaining four Monte 
Carlo integrations. It might also be easier to perform the mass shift calculations in 
position space as only an integration over the relative distance r would be required. In 
momentum space, one has an additional integral as each wave function depends on a 
different three-momentum p. 

In this paper we derive the static QCD-potential in position space with quark 
masses through two loops. We use the known results in momentum space and 
derive the coordinate results through a Fourier transformation of reconstructed ap- 
proximate analytical momentum space expressions. The latter step is necessary due 
to the complexity of the results of Ref. |§ . 

We begin, however, by recalling the definition of the potential through the man- 
ifestly gauge invariant vacuum expectation value of the Wilson loop. Fig. [j] displays 
the Wilson loop W-p = (0|TrPexp {ig § T dx^A^T a ) |0) of spatial extension r and large 
temporal extension T with gluon exchanges indicated. The path-ordering is neces- 
sary due to the non-commutativity of the SU(3) generators T a . In the perturbative 
analysis through two loops considered here, all spatial components of the gauge fields 
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Figure 1: The Wilson- loop T with large temporal extent (T — > oo) from which the 
position space potential is defined. Through two loops in four dimensions, gluons 
connecting the spatial source lines can be neglected. 



A£(r, ±T/2) can at most depend on a power of logT and are thus negligible here. 
Furthermore, Wr ' 1 —^> exp(—iTE ), where the ground state energy E is identified 
with the potential V. Thus we arrive at the definition: 

V(r,m) = - lim — log(0|TrP exp (ig I dx.A^T 1 ) 10) (1) 
t^oo %T \ Jr J 

Writing the source term of the heavy charges, separated at the distance r = |r — r'|, as 

j;(x)= W T a [5(x - r) - 5(x - r')] (2) 

and neglecting contributions connecting the spatial components, the perturbative po- 
tential is given by 

V(r,m) = - hm log(0|TrT exp (j d'xA^J^x^j |0) (3) 

In the above equation t> M = 5^o due to the purely timelike nature of the sources. For the 
same reason the path-ordering is replaced by the time-ordering operator T. Expanding 
Eq. H] perturbatively we find the position space Feynman rules for the source-gluon 
vertex and the source propagator respectively: igT a v^ and —i9(x — £q)5(x — x'). 



The potential in momentum space is the Fourier transform of V(r). It can be calcu- 
lated directly in momentum space from the on-shell heavy quark anti-quark scattering 
amplitude (divided by i) at the physical momentum transfer q, projected onto the 
color singlet sector. The momentum space Feynman rule for the source propagator is 
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identical to the one in HQET |K|: v .^ +ie - For anti-sources, v — > —v is prescribed (or 
the time-ordering reversed). Fig. ^summarizes the momentum space Feynman rules. 
In analogy to HQET, double lines denote the heavy source terms. 

The potential can be used to define the effective charge ay ( the so-called V-scheme) 
through: 



V{Q,m) = -4tcC f 



a v (Q,m) 
Q 2 



V(r, m) = — Cj 



a v (r,m) 



(4) 



where 



-q 2 and both expressions above are related through a Fourier- 



transform. The results for QCD-corrections including massless quarks have been cal- 
culated in Ref . [ 11 1 , and in Ref . jl2] an independent approach found a disagreement in 
the pure glue part of the original results. As both authors agree now on the correctness 



of 1 12], we use in momentum space: 
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where /3 



34, 



C| - fC A T F n f - ACpTprif and C A = 3, C F 



and Tp = | in QCD. The number of massless flavors is denoted by rif. The /5-function 
is here defined as 



— L^rn 



a s (jj?) dlog/i 



n=0 



47T 



For the case of massless quarks, the first two coefficients (3q and (3i are renormaliza- 
tion scheme invariant and the /3-function is gauge invariant to all orders in minimally 



subtracted schemes [13 . 



In coordinate space we have in the massless case fll4] 
a v (r) 



7T 



7T 



(9) 



v,k 



i/v.k 8 a , b 



igTa,b v T 



Figure 2: The momentum space Feynman rules used in the calculation of Ref. 0. The 
^-prescription is analogous to the conventional quark propagator. For anti-sources 
v — > —v must be used. 



with 



r>i(r,/x) 
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where r' = rexp(7s). 

From the renormalon point of view, possible power corrections in momentum space 
can at most be of the formQ: 

A 2 



V(Q) 



-AttC, 



Q 2 



const, x ®R D + ... 



(12) 



which follows from Lorentz invariance (since v.q = 0) and dimensional arguments |16 



In this connection it is interesting to note that there is a close connection between the 
potential and the pole mass. Both are affected by the same renormalon ambiguity and 
writing [TjJ: 

"^' )E " + 5jL,p v («l (13) 



1 We still assume that Q > Aqcd and are concerned only with the form of the leading power 
corrections from a renormalon analysis. One cannot learn anything about non-perturbative (i.e. 
confining) effects through this approach. 
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the so-called potential subtracted mass mps(^f), which depends linearly on the cutoff 
[if, can be used as a less IR-sensitive mass parameter for threshold expansions than 
the pole mass m JT7 . 



In position space, however, 

V(r) = -C F ^^- (l + ... + const, x A QCD r + const' x A 2 QCD r 2 + ...) (14) 

and it is thus to be expected that the long-distance contributions to the coordinate 
space potential are parametrically larger than for the momentum space potential. 
Thus, the position space potential is likely to be more slowly convergent. This feature 
could be interesting when comparing the full lattice results with the perturbative so- 
lution at intermediate distances |T5|| . The form of the linear term in Eq. |TJ] is also 
motivating the consideration of the derivative of V(r), i.e. the definition of the strong 
coupling from the force between the two heavy quarks [0. 

Before considering the effects of massive quarks in the quantum corrections to the 
potential, a few general remarks are in order concerning the higher order perturbative 
behavior. The effective theory used for the calculation amounts to a non-local approach 
with Q ~ Mv , where M denotes the heavy quark mass, i.e. where the gluons are 
always kept off-shell. Through power counting arguments one can see that through 
two loops only those gluons need to be considered. At the three loop level, however, 
on-shell gluon contributions of order Q ~ Mv 2 cannot be omitted and a treatment 
according to the standard definition would lead to an infrared divergent potential |19 



These ultra soft terms are expected to be canceled, however, by additional diagrams at 
that order. In particular, new diagrams connecting also the spatial components of the 
Wilson loop in Fig. p] will contribute at three loops. In Ref. [^0] it is indicated how the 



problem already shows up at the two loop level in three dimensions and in Ref. JT9 
an infrared save definition at higher orders is suggested. In four dimensions, however, 
and to the order we are working here no such problems occur and the definition |3| is 
infrared safe. 

The paper is organized as follows. In section |2] we review the results of the mass 
corrections to the heavy quark potential in momentum space at the two loop order. 
From the physical Gell-Mann Low equation we reconstruct a simple analytical approx- 
imate expression for the one- and two-loop coefficients. These results are then used to 
obtain the position space potential in section |^. In section [| we then discuss the effect 
of massive quarks on the force between two static sources and close with concluding 
remarks in section [|. Explicit results for the coupling definition through the force in 
coordinate space are given in the appendix, where we also discuss briefly the effect of 
mixed massless and massive quark loops at the two loop level. 
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2 Momentum space results 



The Monte Carlo results of Ref. || can be used to obtain the two loop scale de- 
pendence of the static QCD potential in momentum space |J. The difference to the 
conventional /3-function in the Callan-Symanzik pT| , p2[ approach is that the physical 



quark anti-quark system is governed by the exchanged momentum and independent 
of the renormalization scale to each given order. Thus we follow the Gell-Mann Low 
approach and using a running MS-mass m(fi) we have: 



a v (Q,m(fjL)) = a m {v) ( 1 +v 1 (Q,m(ji), /jl) 



7T 



7T 



1(15) 



where the massless limit of the coefficients V\ and t>2 is given in Eqs. [] and [7]. The 
physical charge av(Q,m(fj,)) cannot depend explicitly on the renormalization scale \i 
and the explicit /i-dependence on the right-hand side of Eq. (|15|) cancels to the order we 
are working. Fig. |^ gives the Feynman diagrams for the fermionic contributions to the 
two-loop coefficient V2{Q, m(fi), //). The mass-counterterm chosen for the Feynman 
diagram labeled gse^ determines the mass-parameter which has to be used in the 
one-loop coefficient V\(Q , m(n) , //) . In Ref. || we considered the flavor-threshold 
dependence of heavy quarks and related the running mass to the pole mass which is 
renormalization-scale independent and gives explicit decoupling. This also provides a 
physical picture as well as a straightforward Abelian limit. 



The relation between the MS mass and the pole mass m is given by |29 

m(/i) = m 
Inserting Eq. ( |16D into Eq. (|15|) gives 
a v (Q,m) = «Ms(/i) 



1 - C F 1 + - log — 

7r V 2 m 



1 + Vi{Q,m,jJL) h 



(16) 



7T 



a— (a) 

h(Q,m,/i) + A m (Q,m,/i)]^4^ 



7T Z 



(17) 



where A m (Q,m, fx) denotes the contribution arising from v\ when changing from the 
running mass to the pole mass: V\(Q , m{ji) , fi) = t>i(Q, m, ji) + A m (Q, m, /i) atA ^ . 

The Gell-Mann Low function |23|] for the ^-scheme is defined as the total logarith- 
mic derivative of the effective charge with respect to the physical momentum transfer 
scale Q: 

(Q\_ da v (Q,m) _^ M o# 2 {Q,m) 
* y \m) = d\o g Q = h ^ ^ ' (18) 
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gse,: gse 2 : gse 3 : 




Figure 3: The two-loop massive fermionic corrections to the heavy quark potential in 
the Feynman-gauge (from |§). Double lines denote the heavy quarks, single lines the 
"light" quarks with mass m. The first two rows contain diagrams with a typical non- 
Abelian topology The middle line includes the infrared divergent "Abelian" Feynman 
diagrams. They contribute to the potential only in the non-Abelian theory due to 
color factors oc CpCa- In addition, although each diagram is infrared divergent, their 
sum is infrared finite. The infrared finite Feynman diagrams with an Abelian topology 
plus the diagrams consisting of one-loop insertions with non-Abelian terms are shown 
in the last two rows. 
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where in the massless case the coefficients ijffl and ipy are given by, 

^ 0) (m = 0) = ^C A - \t f N f = H _ 1 Nf , (19) 
o 6 2 6 

^(m = 0) = ^Ci - \C A T F N F - \C F T F N F = h l - ^N F . (20) 

For the massive case all the mass effects are absorbed into a mass-dependent function 
Np. In other words we write 

- t4<0 w 

_ 51 19 (Q 



\m/ 4 12 ' \m/ 

where the subscript V indicates the scheme dependence of N FV and N FV . 

Taking the derivative of Eq. fllTf) with respect to log Q and re-expanding the result 
in ay (Q,m) gives the following equations for the first two coefficients of \EV: 

( 0) /Q\ = _d» 1 (Q,m,ri 
\mJ dlogQ 

ih m (Q) = rf N(Q,m,/x) + A m (Q,m,n)] | m d Vl (Q,m,fx) 

v \mj dlogQ ' dlogQ 



Eqs. (|23|) and (pl|). Rather, and are functions of the ratio of the physical 



momentum transfer Q = \/—q 2 and the pole mass m only The expression for if) 



The argument Q/m indicates that there is no renormalization-scale dependence in 

and 

(0) 
V 

agrees with our result in Ref. [21]. In Eq. ( [24] ) the derivative of the A m (Q, m, /i)-term 
comes from using the pole-mass instead of the MS mass whereas the remaining mass 
dependence in Eq. (53) is arbitrary in the sense that a different mass scheme is formally 
of higher order. In addition we note that the contribution 2v\dvi/ dlogQ cancels the 
reducible contribution (labeled 2vp in Fig. ^) to v-i\ it is thus sufficient to consider one 
quark flavor at a time for the two loop Gell-Mann Low function. In the appendix we 
describe how to treat the effect of massless quark loops (u,d and s) in mixed Feynman 
diagrams on the amplitude level as well as the case for two different mass flavor loops. 

Because of the complexity of the integrals encountered in the evaluation |§ of the 
massive two-loop corrections to the heavy quark potential, the results were obtained 



numerically using the adaptive Monte Carlo integration program VEGAS [25|. Thus 



the derivative of the two-loop term v% was calculated numerically, whereas the other 
terms in Eqs. (23) and (^4j) were obtained analytically. The results are given in terms 
of the contribution to the effective number of flavors Np V and Np V (jjj^J in the 
V^-scheme from a given quark with mass m defined according to Eqs. (pTT) and (p2|) 



respectively. The Appelquist-Carazzone |26| theorem requires the decoupling of heavy 
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masses at small momentum transfer for physical observables. Thus Npy (&) goes to 

zero for Q/m — > 0. The massless result Npy — > 1 is also be recovered for large scales. 

The calculation presented in Ref. |§ required the evaluation of four-dimensional 
scalar integrals. The results in Ref. are based on 50 iterations of the integration 
grid each comprising 10 7 evaluations of the function which where needed to achieve 
adequate convergence. Even so, the Monte Carlo results still are not completely stable 
for small values of Q/m, especially in the light of the numerical differentiation required 
in Eq. Nevertheless, accurate results can be obtained by fitting the numerical 

calculation to a suitable analytic function. 

The one-loop contribution to the effective number of flavors Np follows from the 



standard formula for QED vacuum polarization. In Ref. |24j we used the simple rep- 
resentation in terms of a rational polynomial []30| : 



nP v (2)k = (25) 

Fy \m) 1 + 5.2 f l + a f 

which displays decoupling for small scales and the correct massless limit at large scales. 
Similarly, the numerical results for the two-loop contribution can be fit to the form 

Q 2 Q A 



m) Q j_ Q 



The parameter values eij and the errors obtained from the fit to the numerical calcu- 
lation in the V^-scheme for QCD and QED are given in Ref. |J. Similar decoupling 
forms have been used for interpolating the flavor dependence of the effective coupling 
in the momentum subtraction schemes (MOM) 32| . 



In the case of QCD we obtain the following approximate form for the effective 
number of flavors for a given quark with pole-mass m: 

(-0.571 + 0.221^- 9 

N m (Q\ „ V m / - (27) 

1 + 1.326-^7 + 0.221-^r 




m 2 m 4 



and for QED 



(1.069 + 0.0133^)^ 

< u) w - — & — V • (28) 

V J 1 + 0.402^- + 0.0133-^ 



m? m 4 



The results of our numerical calculation of Npy in the ^/-scheme for QCD and QED 
are shown in Fig. f|. The decoupling of heavy quarks becomes manifest at small Q/m, 
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and the massless limit is attained for large Q/m. The QCD form actually becomes 
negative at moderate values of Q/m, a novel feature of the anti-screening non-Abelian 
contributions. This property is also present in the (gauge dependent) MOM results. In 
contrast, in Abelian QED the two- loop contribution to the effective number of flavors 
becomes larger than 1 at intermediate values of Q/m. We also display the one- loop 
contribution Np V which monotonically interpolates between the decoupling and 
massless limits. The solid curves displayed in Fig. f| show that the parameterizations 
of Eq. ( JZTD which we used for fitting the numerical results are quite accurate. 




-0.5 



10 



10 



Q/m 



(1) 



Fi gure 4: The numerical results for the gauge- invariant Npy 
and QCD (triangles) with the best \ 2 fits °f Eqs. fl28|) and 



in QED (open circles) 
superimposed respec- 
tively (from Ref. ||). The dashed line shows the one- loop Np V function of Eq. (fffi). 
For comparison we also show the gauge dependent two-loop result obtained in MOM 
schemes (dash-dot) |32|. At large £ the theory becomes effectively massless, and 



both schemes agree as expected. The figure also illustrates the decoupling of heavy 
quarks at small — . 

In Ref. |3j it was shown that the Abelian limit displayed in Fig. ^ agrees with the 



well known literature results [|33], |34], |35j and that the full QCD result is independent of 
the renormalization scale. The very good agreement of the exact two loop calculation 
with the relatively simple fitting function for Npy of Eq. |27| makes it possible to 
reconstruct an analytical approximate function for the full mass dependent two loop 
coefficient in the next section. 
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2.1 Reconstructing the momentum space potential 

Starting from the general expression for the Gell-Mann Low function in Eq. [18| we 
can obtain av{Q,m) through integration over logQ. Our goal is to reconstruct an 
analytical function for av(Q, m) in terms of OmsQj,) based on the fitting parameters ao 
and ax, a 2 , a 3 from the approximate one- and two- loop solutions respectively. It should 
be clear that by analytical we mean an expression of known functions depending on 
{Q 2 , m 2 , fi 2 , a , di, a 2 , a 3 } in the spacelike (Euclidean) region. It is understood that for 
a continuation into the timelike regime only the full function calculated in Ref. |9j can 
be used, not the approximate form we derive below. With this in mind we can write 



a v {Q,m) - aMs(fi) 



71 
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1$ d log Q + C (0) 



\fj { y ] - 2v 1 ^ (0) ] d logQ + C (1) 
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The integration constants can be functions of m and \i and are fixed by requiring 
that the correct massless limit is obtained. We find for the corrections with one massive 
quark with pole-mass m: 



vi(Q,m,n) 
v 2 (Q,m,fi) 



-^--C A log- + - 

Cl /4343 2 7T 4 
+ 4tc z 



16 V 162 
+T F 




C A /1798 56 \ C F /55 
"16 V~81~ + ~ 16 \ J 



16Cs ) + 2 ^T F 



+ ■ 



19 



6 



1 . m 4 + a 3 m 2 Q 2 + a 2 Q 4 

A 4 + 

1 a 2 /z 4 



4 



ai 
2 



a? — 4a 2 



:X 



log 



2a 2 Q 2 + 



4a 2 + a 3 )m 2 



2a 2 Q 2 



4a 5 



55^ , Q 2 
+ — C A log + 



a 3 Jm z 



54 



31 _ . Q 2 + a Q m 2 

— C A log - 

54 /r 



Q 2 



10 Q 2 + a m 2 11 2 

— T F log — C A log 

27 /i 2 36 /i 2 



-T^C^log^ lo 



Q 2 Q 2 + a m 2 1 2 Q 2 + a m 2 



18 



+ -T F loj 



(30) 



In terms of the running mass in the MS-scheme the solution reads: 



1 4 + 3 log R 

V2(Q,m(fj,),fj) = v 2 (Q,m,ii) + —C F T F 

6 1 + 



(31) 
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100 150 200 250 300 350 400 



Figure 5: The comparison between the exact results from Ref. || (open symbols) and 
the reconstructed solution in Eq. |31] (solid lines) for the bottom (triangles) and charm 
(circles) quarks. The absolute value of the fermionic contributions (proportional to 
Tp) times is shown. The scale \i was chosen to coincide with the quark masses in 
the MS-scheme of Ref. [[J. It can clearly be seen that Eq. [H] is in good agreement 
with the full analytical result over all perturbative values of the momentum transfer 
Q within the statistical Monte Carlo and fitting errors of a few percent in each case. 



The results are written in such a form that the limit m — ► is obvious and can be 
seen to agree with Eq. |5|. Eq. can be directly compared to the exact calculation in 



Ref. IBfl. The latter was obtained in the MS-scheme, so we must use 




Fig. H shows the good agreement of approximate solution in Eq. ^T] with the full 
result for different input parameters. The figure also displays the fact that the fitting 
parameters ao, ■■, as of Ref. || are optimized for the flavor threshold region Q ~ 0(jn). 
In Fig. |]the two loop Q-dependence of the physical ay(Q, m) charge is compared with 
the massless calculation. In the massive result the bottom quark mass m;, = 4.5GeV is 
used and it can be seen that the two massive results (pole and running mass) differ by 
several percent from the massless (rif = 1) calculation in the region Q ~ m . It should 
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a v (Q,m) 

0.38 | 




Q 

Figure 6: The momentum dependence of av(Q,m). We use \i = M z , ay^{M z ) = 0.12 
and rif = 5 throughout. The solid line is the massless result, while the two dashed 
lines on top of each other are massive results for the b-quark with the pole and MS 
mass. 
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also be stressed that we display the full Q-behavior of ay(Q,m) including constants, 
not just the running according to the Gell-Mann Low function. 



3 Coordinate space results 



In this section we will present results for the coordinate space potential based on the 
Fourier transform of Eq. [JTI It should be emphasized that the Fourier transform of 
av(Q, m) in the strictly perturbative sense does not exist (Landau-pole) and that only 
the expanded coefficients can be used. This will be shown below. In general, we have 
the following relations: 
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(33) 



At fixed orders, there is no Landau-pole in Eq. [33] which can best be seen by inserting 
the expansion in Eq. into Eq. [J^: 
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(34) 



It can easily be recognized that the first term in the series is just the Coulomb potential 
Vq = —Cp 01 ™^ . The loop corrections are then given by writing 
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-2 Ei ( 1, \j a 3 — \j a 2 — 4a 2 m r 



+ 2 Ei 1, \ a 3 + Jai-4a 2 mr 



55 

— C A (log(/x r) + 7 £ ) 



+ 



54^ _ 27 Tf 



log — — + 2Ei (1, ^/oq m r) 



11 

36" 



7T 



112 T 2 

+4(log(/x r) +7i?) 2 ) - — CU-Zi(r,m,//) + -f -J 2 (r, m, /i) 

7 10 7T 9 7T 



where the integral representations are defined as follows: 

roc , 

Ei(l,x) = / exp(-t)- 

J X 



dt 
T 



Xi(r,m,/i) = / log 
Jo 

roc 

l 2 (r,m,/j,) = / log 
Jo 



Q 2 + a m 2 Q 2 . dQ 

~ log — sm(Q r) — 

Qf-i /r Q 

i Q 2 + a m 2 . dQ 
sm{Q r) — 



Q 



(37) 

(38) 
(39) 
(40) 



Eqs. |38], [39] and ^0] are readily calculable numerically, for instance with the mathemat- 
ical package MATHEMATICAL Again, if one uses the MS-running mass m(/x) instead 
of the above used pole mass m the result must be modified in the following way: 



1 



v 2 (r, m(fi),fi) = v 2 (r, m, y) + ^C F T F I 4 + 3 log — I (1 - exp(- v /a^ m r)) (41) 



6 



77T 



The effect of the two loop mass effects in position space are depicted in Fig. [7| for 
the bottom quark with m& = 4.5GeV. The pole- and MS-mass corrections are almost 
identical on the scale (dashed lines) and the effect grows with larger distances. Since in 

2 For Eqs. [39] and |40| the integration method should be adapted to the oscillatory behavior of the 
integrands. 
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dv(r,m) 




Figure 7: The distance dependence of oty{r, m) for the same parameters as in Fig. 
r is in GeV" 1 and 0.2fm « lGeV -1 . 



D.9 - 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 

r 



Figure 8: The distance dependence of ay(r,m) for different choices of the renormaliza- 
tion scale /x. The solid line corresponds to \i = M z as in Fig. [7|. The dash-dotted line 
to /i = - and the dashed line to the "natural choice" /i = — The source separation 
r is in GeV -1 . 
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V ( r > m ) G eV 




Figure 9: The distance dependence of V(r,m) with the same choices of the renormal- 
ization scale /i as in Fig. M. The potential is given in units of GeV, r in GeV -1 . 



the renormalization group (RG) logarithms always enter as log(/xr) + 7_e, the "natural" 
choice is given by jx = In Ref. |27]] it was shown that this scale almost identically 



reproduces the BLM-results in the massless case, thus consistently reabsorbing the 
large RG-group effects. As in this case we also have to evaluate the two-loop running 
of the MS-coupling using: 



An 



v 2 — 



V 



log log 



l QCD 



Pi 



log 



^QCD 



(42) 



/ 



where we normalize the QCD-scale parameter Aqcd such that os^g(Mz) = 0.12 which 
corresponds to Aqcd = 0.25GeV and we keep rif = 5 fixed. Fig. ^| displays the effect 
of the various scale choices on cny(r, m). Although the physical coupling «y(r, m) has 
no renormalization scale dependence to the order we are working, the difference seen 
in Fig. |8| is due to uncanceled higher order terms and, in general, increases at larger 
distances. Substantial deviations can be seen for ay(r, m) > 0.2. The potential is 
shown in Fig. |^ and one has to use the standard conversion factor lfm" 1 w 0.2GeV 
in order to convert the distance to fm. Again, the scale dependence is due to higher 
order contributions. The mass dependence of ay(r,m) is displayed in Fig. pi We 



17 



a v (r,m) 




0.5 1 1.5 2 2.5 3 3.5 4 4.5 



Figure 10: The dependence of av(r, m) on the MS mass parameter m(/z). The figure 
keeps r = 0.1 fm fixed. The diagram is displayed for the "natural choice" of the 
renormalization scale \i = — 1^-. 
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keep r fixed at 0.1 fm and vary the MS-mass parameter between and 5GeV. While 
perturbation theory might not be applicable for distances larger than O.lfm, the effect 
is significant and can amount to a few percent for MS-masses up to lGeV for the 
physical charge a>y(r,m) compared to the massless result. For the bottom quark mass 
the effect is of order 10 % at r = O.lfm. 



4 The force between two static sources 

In this section we discuss the concept of defining the strong coupling from the force 
between two static color-singlet sources 0. In coordinate space the force is simply 
given by: 

„, , dV(r,m) _ a F (r,m) , . 

F(r, m) = ^-^ = -C F FK 1 ' 43 

or r 2 

For large distances there is no required sign change for a F (r, m) and its accompanying 
/3-function is unique. For the massless coupling we can simply write: 

a F (r) = a v (r) (1 - 20 v (a v (r))) (44) 

where j3y denotes the massless (3-i unction in the V-scheme. From Eq. [H] it follows 
directly that the massless relation between the a^-charge and the MS-coupling to two 
loops is given by 

/ \ / A l , , , OmsV) , , ( "msK) /.-v 

a F (r) = a m (r ) I 1 + — - — + f 2 I — - — I (45) 

where r' = rexp(7^) and the constants are given by 
35 1 

fx = ~-C A + -T F n f (46) 



f2 



35 . 1 




-C A + -2 


36 9 


1 


V 7513 


16 


LI 162 




'31 


-( 





^ 2 

27 4 



In general, the relation for the massive case is given by 

2 <9(ay(r, m)/r) 



a F {r,m) = —r 



dr 

2^ 



«nff(A*)|l + /i(r,m I /*)^^ + / 8 (r I r^/*)(^^) ) (48 

at the two loop level. Explicit expressions for f%(r, m, //) and f2{r,m,fi) are given 
in the appendix. The distance dependence of the perturbative coupling definition 
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a F (r,m) 




F ( r . m )ceV 2 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 

r 



Figure 11: The position dependence of ap(r, m) and the force F for different choices 
of the renormalization scale \i. The scale choice corresponding to the respective lines 
is the same as in Fig. 13. The force in the lower plot is in GeV 2 , r in GeV -1 . 
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of ap(r, m) is presented in Fig. |TT] for the natural scale choice (dashed) and fixed 



/i = Mz (solid). The former leads to a lower value of Op(r, m) at larger distances. 
Overall, it can be seen that the coupling defined from the force is much smaller at a 
given source-separation r than azv(r,m) (and even smaller than (^(/i = — ^)). Even 
at r = 0.2 fm (= lGeV -1 ) still seems to be acceptable from a perturbative point of 
view. This could make a comparison with lattice results very interesting. Overall, the 
renormalization scale dependence is also strong for ap{r, m). The lower graph displays 
the force itself in units of [GeV 2 ]. The mass dependence of ap(r, m) is given in Fig. p"2] 
analogously to Fig. [To] by varying the MS-mass parameter between and 5GeV and 
for two fixed values of r. The dependence on the mass terms is not as large for small 
MS in comparison to ay{r, m), however, can still be a few percent. It can also be seen 
that the effect for the MS bottom quark mass is roughly 7% compared to the massless 
result at r = O.lfm, while it is about 10% at r = 0.2fm. 



5 Conclusions 

We have calculated the coordinate space static QCD-potential through two loops with 
quark masses. The result is obtained by reconstructing the exact momentum space 
Monte Carlo results from Ref. [f| by analytically fitting the Gell-Mann Low function. 
The reconstructed results are in good (few percent) agreement with the fermionic 
results of the exact calculation. Based on these solutions we have obtained the Fourier 
transform in coordinate space and the force between the heavy quarks. The Qy(r,m) 
coupling itself grows quicker than a^/i), which is mainly due to its larger value at 
Q = Mz- The mass-dependence is strongest for small MS masses up to lGeV at large 
distances (O.lfm), where the perturbative regime breaks down. At these distances the 
effect of keeping a massive charm or bottom quark can be a significant effect, although 
perturbation theory might not be trustworthy anymore at distances of order 0.1 fm. 
An important indicator of this fact is the strong renormalization scale dependence from 
uncanceled higher order terms which can lead to substantially different results. 

For the corresponding definition of ap(r, m) we find also a strong dependence on 
the mass terms. The effect is not as large for small MS masses, however can still be 
a few percent. At 0.2fm the effect is of order 10% for the bottom quark mass. The 
renormalization scale dependence from uncanceled higher order terms is also large for 
af(r, m). The overall smallness of this coupling parameter, however, indicates that it 
could serve as a suitable quantity to compare the perturbative and non-perturbative 
approaches even above scales of (9(0.1 fm). 

In terms of practical applications, the relatively simple form of Eqs. |30] and j35[ 
makes it possible to investigate the uncertainty of a massive charm contribution in the 



determination of the bottom quark mass [36 . 
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0.94 




0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 

m 



Figure 12: The dependence of ap(r, m) on the MS mass parameter m. The figure on 
the top keeps r = 0.1 fm fixed while the one on the bottom has r = 0.2 fm. Both 
diagrams are displayed to the "natural choice" of the renormalization scale \i = — 7^-. 
The effect of mass terms is clearly visible. 
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Appendix 

Explicit two loop results with quark masses for ap(r, m) 

In this appendix we give the explicit expressions for the two loop results for the strong 
coupling parameter defined from the force between two static color singlet sources. 
The general definition is given in Eq. 48 and is given by: 



, 35 11 

fi(r, m, fi) = C A [ -— + — (log(/i r) + 7s) ) + T F 



5 2 , . — . 
-- + -exp {-y/ao m r) 



+ 



log— — h 2 Ei (1, t/oo m r) 



(49) 



r<2 
16 



7513 229 



484 /, / \ X 

+— (log(^ r) + 7^) 



162 27 

2^ 



7T 4 22 

f 2 (r, m, n) = ^ ( + - — + —( 3 - ^ (log(/x r) + lE ) 
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Figure 13: The Feynman diagram which can contain both massless and massive (thick 
line) quark loops or two flavors of different non-zero mass at the two loop level. 



-2Ei [ 1, \/ a 3 + \/al — 4a,2 m r 



55 

— CU(log(/x r) + 7jB ) 



~^C A 4(log(H+7.) +y -YV Jl + 9V X2 



where 



Q2 _|_ aQm 2 q2 

log — cos(Q 



/■OO 

J 3 (r,m,/i) = / log 

M{r,m,fi) = / log 

Jo u z 



cos(<5 r)dQ 



For the running MS-mass one needs to add the following contributions: 
f 2 {r,m(ji),fi) = / 2 (r,m,/i) - -C F T F [4 + 31og^ ) {^/a^m r exp( 



(50) 



(51) 
(52) 



6 

1 -exp(-v^mr))) 



do m r) 



(53) 



Mixed massive and massless fermionic contributions at the two 
loop level 



In this part of the appendix we briefly comment on the type of contribution depicted 
in Fig. [13]. This diagram can contain the same flavor in both loops as well as different 
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ones. In the former case, the contribution is already contained in the discussion above. 
In section [2.1| we demonstrated that for the Gell-Mann Low function only one flavor 
needs to be considered. On the level of the potential, however, we have also the latter 
case represented in Fig. [13. 

In a practical situation, we can often set the masses of the three lightest quarks, 
u,d and s to zero. In this case the Feynman graph [13| contributes with a multiplicity 
factor two and we find the resulting expression: 

4vrCV (31^ 5 m 1 „ , u 2 \ 

2 x F,g.l V, = -2-^- [-C A - -Tpnf + jA log iL j x 

[36 Ca ~ 9 Tf + ~2° A bg Q 2 - ^ QT^^ ) — (W) 

In the case of two non-zero and different quark loops we find analogously: 

4vrC F /31^ 5 m 11 , u 2 1 , u 2 \ 

2xFig.l3 m , ro , = -2— — C4--TP + — C A log^--Tplo f ] - 



g 2 V 36 9 12 3 Q 2 + a m 2 J 

(31 5 11 /i 2 1 /x 2 \ aUv) ,^ 

In the above expressions we used again the approximate one loop solution for the 
vacuum approximation. If desired, it can of course be substituted with the exact 
function containing more complicated expressions ||24)| . We also already divided by a 
factor i to obtain the contribution to the potential. 
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